Quantum wells, wires and dots with finite barrier: 
analytical expressions of the bound states 
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From a careful study of the transcendental equations veri- 
fied by the bound states energies of a free particle in a quan- 
tum well, cylindrical wire or spherical dot with finite potential 
barrier, we have derived analytical expressions of these ener- 
gies which reproduce impressively well the numerical solutions 
of the corresponding transcendental equations for all confine- 
ment sizes and potential barriers, without any adjustable pa- 
rameter. These expressions depend on a unique dimensionless 
parameter which contains the barrier height and the sphere, 
wire or well radius. 

PACS numbers: 73.21.Fg 



The present experimental and theoretical efforts in 
semiconductor physics are essentially devoted to confined 
structures: quantum wells, wires or dots. In most calcu- 
lations, the potential barrier between the confined struc- 
ture and the outside semiconductor is assumed to be in- 
finitly high for simplicity. When the barrier is finite, the 
number of bound states, which is infinite for infinite bar- 
rier, can be reduced to 1 or even for spherical dots, each 
energy level being substantially modified by the usual ex- 
perimental confinements. 

The great advantage of the infinite barrier assumption 
lies in the fact that the energies and wave functions have 
analytical expressions so that all theoretical calculations 
using them are rather easy. The introduction of the finite 
barrier height Q-j^] makes everything much more compli- 
cated. Even in the simplest case of one confined direction 
only, namely the quantum well, the bound states ener- 
gies are given through a transcendental equation. Its 
solutions can be obtained either by a numerical calcula- 
tion or by " reading" their values on a set of curves which 
makes their use quite unconvenient in quantitative calcu- 
lations. The cylindrical wire and spherical dot confine- 
ments being obviously more complicated, the situation is 
even worse for the energy of these 2D and 3D geometries. 

The purpose of this communication is to provide an- 
alytical expressions for the energies of any bound state 
level of these three geometries, valid for all confinements , 
i.e. all barrier heights and confinement sizes. We will 
show that these energies have the same analytical form, 
which depends on a unique dimensionless confinement 
parameter v, in which enter the barrier height and the 
confinement extension. Let us add that, once the energy 
levels are known, it is straightforward to get from them 



the exact wave functions of all these bound states. 
1. Results: 

We consider a particle of mass m confined in a sphere 
or cylinder of radius R or in a quantum well of width 
2R, the energy barrier being V. From R and V, we can 
construct the dimensionless parameter v which rules all 
the physics of these finite barrier problems, namely 



V = h 2 v 2 /2mR 2 



(1) 



The infinite barrier limit corresponds to v infinite. From 
this parameter v, we can already note that, as the physi- 
cal scale for the barrier height is ft 2 /2mR 2 , a given barrier 
V between two semiconductors can appear as high or low 
depending on the confinement extension: the larger the 
R, the better the V = oo approximation for the same V. 
In other words, finite barriers effects are going to be very 
important for strongly confined systems. 

Following Eq.(fJ), we are led to measure the particle 
energies in the same unit as V , namely 



E = h 2 a 2 /2mR z = V - ti z z /2mR 



2o2 



(2) 



Note that for ft ~ 0, the energy is close to the top of the 
well while for ft ^ v, a is close to so that the level is 
deep inside the well. 

From the transcendental equation verified by the var- 
ious energy levels of these confined geometries given be- 
low, we find that the number of bound states is controlled 
by the position of the parameter v with respect to a set of 
values Vmin at which a bound state gets out of the well, 
i.e. disappears. A careful study of this transcendental 
equation shows that the parameter a of the bound state 
levels behaves as a max /(l + v^ 1 ) when v — > oo, and as 
v m i<n+"({v— v-min) when v — > p m i n - The three parameters 
{vraim ot max , 7) of these asymptotic behaviours depend 
on the geometry of the confinement and the level under 
consideration. 

• For quantum wells, the levels are characterized 
by one quantum number n, with n — 1,2,3, the 
corresponding parameters being v m i n — (n — l)7r/2, 
u m ax = rim I '2 and 7=1. 

• For cylindrical wires, the levels are characterized by 
two quantum numbers n and ±m with m = 0, 1, 2, • • •. 
For the (n, ±m) levels, the parameter 7 is equal to 
l/sup(m, 1), the parameter a max is equal to the n th zero 

z~m of the Bessel function J m and the threshold v m in for 
level disappearance is equal to z^-i f° r m 7^ 0, z[ n ^ for 
(rn = 0,n^ 1) and for the ground state (m = 0, n = 1). 
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• For spherical dots, the levels are characterized by 
three quantum numbers (n, Z,±m) with I = 0, 1,2, 
they are degenerate with respect to to. For I = 0, the 
parameters are v m in = (n — l/2)n, a max — nix and 7 = 
1, while for I ^ 0, they are is mm = 2j_i/2> a max = 4+\/2 
and 7= 1/(Z + 1/2). 

From the numerical resolution of the transcendental 
equations verified by the energy levels of these three con- 
fined geometries, we have checked that the energy param- 
eter a is amazingly well reproduced, for all levels and all 
confinement is, by the same function 



,-l)+(y-i/ min )(l+(7-l)a mQ . ai /i/ min ) 

(3) 



which is constructed to give the two first terms of the 
a behavior for both v — > 00 and is — > is m i n (see Fig. 
|l]||[|) ■ The solid lines which correspond to the numerical 
solutions of the energy transcendental equations are hard 
to distinguish from the dashed lines which correspond to 
Eq.(||): the discrepancy is indeed extremely small. 

In the particular case of the n th level of a quantum 
well (which exits for is > (n— l)7r/2 only), this gives 



aW(i/) 



v + l 



(n/2-iy 



(4) 



(tt/2 - 1) + (is - (n - 1)tt/2) 



while for the n th level of the / = states of a spherical 
dot (which exists for v > (n — l/2)7r only) this gives 



a\J Q (v) ~ T17T - 



v+l 



2. Derivation: 



(n/2-iy 



(5) 



(tt/2 - 1) + (v - (n - 1/2)tt) 



Let us now outline how we have derived the above re- 
sults. The Schrodinger equation of a particle of mass 
m in the confined geometries considered here reads, in 
terms of the reduced variable £ = z/ R for quantum wells, 
£ = (\p\/R, f) for cylindrical wires and £ = (\r\/R, 9, >p) 
for spherical dot, as 



(6) 



with «(£) = a 2 for ||j < 1 and v(£) = ~(3 2 for ||j > 1. 

The bound states of these confined geometries corre- 
spond to (a, f3) real and positive a and (3 being linked 
by: 



(3 2 = v 2 



(7) 



due to Eq. 



a) ID case: quantum well 



In the ID case, the solution of Eq.(g) is elementary: the 
wavefunction which cancels at ±00 and has continuous 
derivatives at ±1 can be written as: 



a 2ia iai 



+ e 



2i arctan - 



xjj(-l < £ < 1) = a 

< -1) = V(-l)e^ ( « +1) 
> 1) = iP(l)e-^-V 

with a and j3 linked by Eq.(Q) and by 

exp 4 i a = exp(4 i arctan (3 /a 



(8) 



(9) 



Eq. (H) is not one of the standard forms of the quantum 
well transcendental equation found in usual textbooks 
G10|. It is however the most convenient one to extract 
the asymptotic behaviour of a as we now show. 

(i) When v — > 00, the finite values of a corre- 
spond to (3 ~ v infinite. For such v, [3/a is large so 
that arctan/3/a ~ 7r/2 — a/(3. Eq.(||) then leads to 
4a - 4(tt/2 - a/u) + 2(n - 1)tt i.e. a ~ 2f7(l + z/" 1 ) 
with n = 1, 2, 3, • • •. 

(ii) The solutions for the bound state disappearance, 
i.e for an energy level close to the top of the well, cor- 
respond to gl — v and /3 ~ 0. [3/v is then small so 
that Eq.(||) gives (3 ~ 1/(1^ — (n — 1)tt/2) i.e. a ~ 
2/ + - (n - 1)tt/2) 2 ) due to Eq.(@). This gives the 
values of a max , v m i n and 7 listed above. 

In order to check the validity of the approximate 
ct^(v) given in Eq.(|j), it is not necessary to solve 
Eq.(n,ffl numerically: indeed, while these equations do 
not give a in terms of v, they do give is in terms of a as 



is = a/ cos(a — (n — l)7r/2) 



(10) 



The corresponding curves as well as the various a'"' (is) 
are shown in Fig. [l| The fit is impressive; the two sets of 
curves are hard to distinguish, the largest discrepancy is 
smaller than 4%. 



b) 2D case: cylindrical wire 
In the 2D case, we have A|- 



lJLpJL _i_ 1 d 2 



| d^di n C 2 dip 2 

that the wave functions write "0(0 = fm{Q e ±tmip with 
to = 0, 1, 2, • • •. The differential equation obtained for f m 
is just the one of Bessel functions (J m ,l^n) for £ inside 
the well and (J m , K m ) for ^ outside. The solutions which 
are finite for £ = 0, tend to for £ — * 00 and have a 
continuous derivative for £ = 1 read: 



f m (H<l)=AJ m (aO 
K m \P) 

with a and /3 linked by Eq.(Q) and by 

(a m ) _ f3 m K m -i(/3 m ) 



(11) 



(12) 
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In order to get the above exp ression, we have used J_i 
-Jl and K-i = Kx and @: 



(a) J' v = J v -\ — vjzJ v 
{h)K' v = -K v _ x -v/zK v 



(13) 
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FIG. 1. The energy parameter a as defined in Eq.(Q) for 
the n th bound states of a quantum well as obtained from the 
exact expression of Eq.(^) (solid line) and from the analytical 
expression Eq.(Q) ( dashed line). The confinement parameter 
v depends on the barrier height V and half width R through 
Eq.fl). 

(i) When v — > oo, the finite values of a m correspond to 
(3 ~ v infinite. The K ratio being 1 in this limit, we find 
J m (a m )/a m J m -i(u m ) ~ -1/v ~ 0. The values of a m 
are thus close to the zeros of J m . For these a m . we find 



using Eq.(|T3|.a) that the inverse of the LHS of Eq.(p^) is 

close to (a m - Zm )/Jm' so that a m ~ Zm . /(I + i^ 1 ). 

(m) The solutions close to the top of the well corre- 
spond to /3 m ~ and a rn ~ ;a For these /3 m ~ the 



RHS of Eq.(jl2h is close to l/ln/3 for m = 0, 1 ln/3 for 
m = 1 and — j1(m — 1) for to > 2. Being close to zero 
in all cases, we conclude that a m must tend to the zeros 
of J m -i. In the particular case of m = 0, we find that in 
addition to the zeros of J_i, i.e. the zeros of J%, cxq can 
also tend to for v — > 0. Since in this /3 ~ i/ limit, a ~ 



(n) 

'm- 
,(») 



-™i_i(a- 



due 



» 



: )/m 
,„,,,) with 



while the LHS of Eq.(|12|) is close to 

to Eq. (|l3|a) , we get for to > 1 , a m ~ z^'_ x + {v—z. 
while for to = 0, we find ao — fmi« + (f — f-, 

"min = or z\ . 

We thus recover the values of (^Vram, QW:r> 7) listed 
above. Figj^ shows the results of the numerical res- 
olution of Eqs.(Rll^) as well as the value of 04' de- 



duced from the corresponding 



, atmaxyj) inserted 



in Eq.(]3j). The fit is excellent for all confinements, the 
two curves being hard to distinguish, with an exception 
for the ground state at intermediate v: the approximate 



a m=o^ i s slightly below the numerical a, the largest dis- 
crepancy being 14%. 




FIG. 2. The energy parameters a for the n th level of 
the ±m states of the cylindrical wire confinement, as given 
by the numerical resolution of Eq.(Q^) (solid line) and by 
the approximate analytical expressions deduced from Eq.(H) 
(dashed line). The horizontal dashed lines correspond to 
4 X) ^ 2.40, ~ 3.83, 4 1 ' - 5.13, z ( 2) ~ 5.52, z[ 2) ~ 7.01 
with z„' being the n th zero of the Bessel function J m . The 
confinement parameter v depends on the barrier height and 
the cylinder radius R through Eq.(hl). Note that two curves 

z i (with different slopes). 



v (") 



. jil=2 and o^Iq start each at v 



c) 3D case : spherical dot 



In the 3D case, we have A^- = |r J|-£ 2 J| — L 2 so 

that the wave functions write ip{£) = fi{£)Y^ m (9,ip), 
with / = 0, 1,2, ••• and —l<m<l. By setting 
fi = ffz/V?, we find that the differential equation for 
gi is just the one of Bessel functions (J;+i/2) ^+1/2) f° r 
£ inside the well, and (-^+1/2,-^+1/2) f° r £ outside. The 
solution which are finite for £ = 0, tend to for £ — * 00 
and have a continuous derivative for £ = 1 read: 



m < 1) 



/;(£>!) = 



A 

71 



Ji+i/2(aQ 



A J J+ i/ 2 (a) 



^1+1/209 



VlK l+1/2 {(3) 
with a and j3 linked by Eq.(Q) and by 



Jl+i/%{oti) 



(14) 



(15) 



for I > 1 (due again to Eq.(jTJ)). For the I = 
states however, it is simpler to use the explicit expres- 
sion Jy,(x) = (JL)V2 sinx and K 1/2 (x) = {^) 1/2 e~ x 
in Eq. (jlj) . The continuity of /q at 5 = 1 then leads to 



ao/tancvo = ~Po 



(16) 
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so that due to Eq.(0), the transcendental equation for the 
/ = levels of a spherical dot is simply 



a sinao 



(17) 



with (n — l/2)ir < a < nir in order to have (3q > 0. 
Let us start with these I = states. 

(i) For v — > oo, Eq.jT^) leads to | sinao| — > so that 
ao ~ n7r. The expansion of the RHS of Eq.(|l7|) close to 
these maximum values of a gives cxq — n7T / (1 + v ~ 1 )- 

(ii) For Pq ~ 0, Eq. (|l6|) leads to cosao — 0; so that 
«o — — 1/2)tt only in order to have /3q > 0. The 
expansion of the LHS of Eq.([l^) close to these values 
gives a ~ f + - (n - 1/2K) 2 ). 

We now turn to the I > 1 states. 

(j) For — > oo and a; finite, i.e. (3i ~ in- 
finite, the ratio of if is equivalent to 1, so that 

^+i/2(a0/ a i^-i/2(«i) — ~ °- a i is thus close 

to a zero of Ji+i/2- Using Eq.(|13|a), the expansion of the 



inverse of the LHS of ( |15| ) leads to (a; 



Jn) 

"1+1/2 



I + I/2) I "1+1/2 



/(I 



so that ai^o : 
(m) For A 

goes to zero as — (3f / (21—1) so that a; is close to a zero of 



0, and ai ~ 2/ finite, the RHS of Eq.© 



Jl-i/2- Close to these values, the LHS of Eq.(15) is equiv- 



alent to Zj_i J2{ a l ~ 

limit we have a;^o 



» 



-1/2 
(n) 

z ;-i/2 



) due to Eq. (|l^a) so that in this 



+ (!/■ 



(n) 

Z l-l/2 



)/(* 



This just corresponds to the values of 



listed above. FigJ 
resolution of Eq. 



1/2). 



y vain ; ^raax 



shows the results of the numerical 
. 15|) or (0) as well as the values of 
a deduced from Eq.(|3|) or given in equation (||). Here 
again, the fit is increadibly good for all confinement and 
all energy levels. 

d) One specific example: 

The above results are of course given in reduced units 
in order to be universal. Let us apply them to one specific 
case: GaAs/Ga x Ali- x As. For usual x's, the electron 
mass is m ~ 0.07 and the energy barrier for electrons is 
of the order of V ~ 300 meV. For quantum well, wire 
or dot of radius R ~ 40A, Eq.ffl) gives the confinement 
parameter v ~ 3.0. Figure (|l]]3|) show that two bound 
states exist for quantum well or wire, but one only for 
quantum dot. Moreover the energies of these levels are 
significantlty lower than the energies of the corresponding 
states for infinite barrier: in the case of a quantum well of 
half width 40A, Eq.(§) gives a ~ 1.1 and 2.24 instead of 
7r/2 and 7r which corresponds to 45 and 170 meV instead 
of 85 and 340 meV. 

In conclusion , we have given analytical expressions 
for the bound states energies of quantum wells, cylindri- 
cal wires, and spherical dots, which reproduce impres- 
sively well their numerical values, for any level, barrier 
height and confinement size, without adjustable param- 
eters. From them, we can easily obtain the exact wave 
functions of these bound states by inserting the corre- 
sponding values of (a, 0) into Eq.(|^JTl|,|l4]). 




FIG. 3. The energy parameters a for the n* level of 
the (l,m) states of a spherical dot as given by the numerical 
resolution of Eq. (|?],[l5|) or (^) (solid line) or by the analytical 
expression deduced from Eq.(efeq3) or given in Eq.Q. The 
horizontal dashed lines correspond to tt/2, tt, ZgV 2 ~ 4.49, 



37T/2, 



,(1) 



5.76, z, 



(i) 



6.99 with z 



1+1/2 



being the n zero 



"5/2 ~7/ 2 

of the Bessel function Jj+i/2- The confinement parameter v 
depends on the barrier height and dot radius through Eq.(|l|). 
Note that for this 3D confinement, the confinement parameter 
has to be larger than 7r/2 for a bound state to exist. 

* Laboratoire associe au Centre National de la 
Recherche Scientifique et aux Universites Paris 6 et 
Paris 7. 
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